HIGH-FREQUENCY AVERAGING IN THE SEMI-CLASSICAL 
SINGULAR HARTREE EQUATION 
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LOUNES MOUZAOUI 



Abstract. We study the asymptotic behavior of the Schrodinger equation 
in the presence of a nonlinearity of Hartree type in the semi-classical regime. 
Our scaling corresponds to a weakly nonlinear regime where the nonlinearity 
affects the leading order amplitude of the solution without altering the rapid 
oscillations. We show the validity of the WKB-analysis when the potential 
in the nonlinearity is singular around the origin. No new resonant wave is 
created in our model, this phenomenon is inhibited due to the nonlinearity. 
The nonlocal nature of this latter leads us to build our result on a high- 
frequency averaging effects. In the proof we make use of the Wiener algebra 
and the space of square-integrable functions. 



1. Introduction 

In this paper we are interested in studying of the Schrodinger equation in the 
presence of a Hartree type nonlinearity 



(1.1) 



iedfu'^ + — Au"^ = £A(A' * 



'■X 



withe> 0, u" e5'(/xR'*,C),de {1,2,3},/ interval of M, A G K, and /-^(a;) = ^ 
with < 7 < d. This model is of physical apphcation. It appears in many 
physical phenomena like in describing superfluids (see [3]). The presence of e in 
this equation is to show the microscopic/macroscopic scale ratio. For small e > 0, 
the scaling of (jl.ip corresponds to the semi- classical regime, i.e. the regime of 
high-frequency solutions u"^ (t, x) . This system contains the Schrodinger-Poisson 
system (d = 3, 7 = 1), a model which is used in studying the quantum transport 
in semi-conductor devices (see IH[251[5S|). Many authors have interested to study 
the asymptotic behavior of (jl.ip in this case like in [ini IHl US] where it is made 
use of Wigner measure techniques. Several papers like [TJ [TU] are devoted to study 
the strong asymptotic limit of the solutions of (|1.1[) with Hartree nonlinearity in 
the case of a single-phase WKB initial data in the form 

wg(a;)=a'^(a;)e*^(^)/^ 

with ip{x) G C £-independent and a^(a;) G M. 

Due to the small parameter e in front of the nonlinearity, we consider a weakly 
nonlinear regime. This means that the nonlinearity does not affect the geometry 
of the propagation (see Section 3), and technically means that it does not show 
up in the eikonal equation, but only in the transport equations determining the 
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2 LOUNES MOUZAOUI 

amplitudes. The object of this paper is to construct an approximate solution uj„p 
for the exact solution m^ of (|l.ip subject to an initial data of WKB type, given by 
a superposition of e-oscillatory plane waves, i.e. 

(1.2) <pp|^^Jx) =^.g(x) =5]a,(x)e-^-/^ 

j<£N 

We begin the paper by showing the existence of an exact solution to (|l.ip in 

L'^iR'^) n WiW^) where W{R'') denotes the Wiener algebra 

WiR") = nL'iR")) = {/ e 5'(R^C), WfWw ■■= WfWmm^) < +00}, 
and where 

In the next step, we construct an approximate solution in the form 

in L^(]R'') n VF(]R'') by determining the amplitudes aj and the phases (f)j, then we 
proceed to study its stability to justify our construction. 
Fourier transform of the potential K is found to be 



K{0 



L^d,7 



(see [3, Proposition 1.29]) so under the general assumptions of [15] where the con- 
sidered kernel K is such that (1 + |^|)if ('^) G L°°(R'^) we can not deduce the well 
posedness of (|l.ip . It is a critical case for our model because we do not know how 
to construct an exact local solution in iy(IR'^) due to the singularity of K around 
the origin (K ^ L°°). 

To define the framework of the amplitudes (aj)jgN we introduce the following 
definitions. 

Definition 1.1. For d G {1, 2, 3} and < 7 < d we define n £ N as follows 

{2 if d^ I or 2 
2 if d = S and 76 [1, 3[ 
3 if d = 3 and -f e]0,l[. 

We define the space 

r(R'') := {f eL^n W, ff'f e l^ n VK, Vt? e W\ \ri\ ^ n], 

equipped with the norm 

ll/lly(K-)= E \\dV\\L^-nw- 

We set 

EiW") = {a = (a,),6N I (aj)jeN G ^'(N,r)}, 
which is a Banach space when it is equipped with the norm 

l|a||_E(Rrf) = 2.^ IWjWy- 
Now, we can state the main theorem of this work. 
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Theorem 1.1. For d G {1, 2, 3} and < 7 < d, consider the Cauchy problem ()1.1|) . 
subject to initial data Uq of the form (|1.2p . with initial amplitudes (aj)jgN £ E{M.'^). 
We assume that there exists 6 > such that 

inf { I Kfc — Km \,k,m £ N,k ^ m} ^ S > 0. 

Then, for all T > there exists C, eo(r) > 0, such that for any e G]0,eo]. ihe 
exact solution to (jl.ip satisfies u' £ i°°([0,T];i^ n W) and can be approximated 
by 

11"" - '"Qppl|L~([o,T];L2niy) < Ce'^ , 
where /3 = minjl, d — 7} and where u' is defined by 

v^Jt,x)=Y,a,(x-tKi)^^'^'^^^'^y^'^-^\'^^^'^^, 

jeN 

with Sj G R defined in p.9p . 

In particular, for d = 3 and 7 = 1 we obtain an approximation result for 
Schrodinger-Poisson equation with /? = 1. This confirnis in a general way what 
was guessed in |^ for the Schrodinger-Poisson system. An important feature of our 
justification consists in the absence of phenomena of phase resonances like we can 
see later. 

2. Existence and uniqueness 
We consider the following Cauchy problem 

(2.1) idtu+-Au = X(K*\u\^)u, 

subject to an initial data uq G L'^{R'^) n W{R'^), with K{x) = t^,0 < 7 < d. For 
< 7 < d, the Fourier transform of K for ^ G M"^ is 



KiO 



^d,^ 



In the following we denote /Ci = iir(^)l[|^|^i] and IC2 = i^(C)l[|5|>i], so, 
K = ICi+IC2 with /Ci e LP(R'^) for aU p G [1, ^[ and IC2 G L9(R'') for aU 
q g]^^, +00]. Let g{u) ^ {K * \u\'^)u. For e > we set 

U%t) = e'"^^, 
with U{t) := U^{t). We recall some important properties of WiW'-). 
Lemma 2.1. Wiener algebra space W{W^) enjoys the following properties (see 

HIS];.- 

i. W{W^) is a Banach space, continuously embedded into L°°{M.'^). 
ii. W{M.'^) is an algebra, in the sense that the mapping {f,g) i— > fg is continuous 
from M^(R'^)^ to VF(]R'^), and moreover 

yf,geW{R''), WfgWw^WfWwMw- 
Hi. For all t G M, U^t) is unitary on WiW^). 
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iv. For all s > i there exists a positive constant C{s, d) such that for all 

\\f\\w^Cis,d)\\f\\H^. 

The following theorem ensures the existence and uniqueness of the solution to 
(EH): 

Theorem 2.1. We consider the above initial value problem (j2.ip with 

uq G L'^{M.'^) n Ty(R'^). Then there exists T > depending on llMoHi^nvy (^"^T-d a 

unique solution u G C{[0,T]; L'^{R'^) f] M^(]R'^)) to (|2T1) . 



Proof. The following lemma will be useful to prove the above theorem : 
Lemma 2.2. Let h£L^r\W,fi,f2&L'^(^W and K like defined as above. Then 

(2.2) \\K*h\\w^\\lC,\\LA\h\\L^ + \\K.2\\LAMw, 

in particular we have 

Il5(/i) - 9{f2)\W-nw ^ Cdl/illi.nw + IIMIi^nwOII/i - MU^nw^, 
for some positive constant C . 

Proof. By definition of the W norm : 

||A-*/i||w = C||(/Ci+/C2)/i|Ui 

<C||/Ci|Ui||/?1|l=»+C||/C2||l-||/J||li 

^C\\lCi\\L^\\h\\L^+C\\}C2\\L^\\h\\w 

Moreover, for fi,f2 & L^ r\W 

\\9ifi)~9if2)\\L^nw = 

= ll(^* (l/lP - 1/21'))/! + (K * |/2n(/l - f2)\\L-nW 

^ m * (l/iP - \f2\'))fi\\L^nw + \\{K * |/2n(/i - f2)\\L^nW 
We control the last two norms. We have 

\\{K^\f2\')ifi-f2)\\L^nw = 

= \\{K*\f2\^Kh - /2)||l^ + ||(A'* |/2n(/l - f2)\\w 
^ ||i^* IM'IIl.^II/i - f2\\L4K*\f2\^w\\fl - f2\\w 

^C\\K*\f2\^\\w\\fi-.f2\\L^nw 

^C||/CiiUl||/2||i.||/l-/2||L^nW + C|1^2|lL^|l/2|Pwll/l-/2||L^nH', 

where we have used (|2.2p and the embedding of T4^(E^) into _L°°(E^). Remark that 

l/iP - I/2I' = ^(/i - f2)iTi + h) + \{Ti-T2){h + h). 



From ()2.2p and Lemma 12711 iii. we have 
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||(i^*(l/il'-|/2p))/i||L^niy< 

^\\KH\fi\'-\f2\')\\L^\\fi\\L^nw + \\K*{\h\'~\f2\')\\w\\fi\\L^nw 
^ C\\K * i\f,f - \hf)\\w\\fi\\L-nw 
^ C\\K * (ifi - f2){Ti +T2mw\\.fi\\L^nw 
+ C\\K*{{Ti-h){fi + f2))\\w\\fi\\L-nw 

^\\)Cl\\L4{fl~f2){Tl+h)\\L4h\\L^nW 

+ C\\IC2\\L^\\{fi-f2){Ti + h)\\w\\h\\L^-nw 

+ C\\IC,\\L4{J-h){h + h)\\L4h\\L-^nw 

+ C\\IC2\\L^\\{Ti-h){.h + f2)\\w\\h\\L^nw 

^ C'll/ilL^nwII/i ^ f2\\L'^nw + C'll/i||L2nH'll/2||L2nwll/i ^ f2\\L'^nW 

+ C||/2|li2nwll/i ~ f2\\L^nw 

^ C'dl/illi^nM/ + ll/slli^nM/)!!/! ^ f2\\L^nw, 
where we have used the Cauchy-Schwartz incquahty. The desired control follows 
easily. D 

Let T > to be specified later. We set 

X^{ue C{[0,T];L^ n W), \\u\\L^^[o,nL-nw) ^ 2\\uo\\L-nw}- 
Duhamel's formulation of (|2.ip reads 

u{t) = U{t)uo -iX U{t - t){K * |up)u(r) dr. 
Jo 
We denote by $(u)(i) the right hand side in the above formula and 
G{u) = $(w)(t) — U{t)uQ. For q > 1 and a space S we define 
hh^^s := \\u\\L^ao,ns) for u e i«([0,T];5). Let ueX. We have 

mu)\\L^L-nL^ ^ mu){t)\\L^L^. + \mu)it)\\L^W: 

and 

^ WuoWl- + C\\K * \uf\\L^L^\\u\\L^L-T. 

We apply Lemma \T^ after replacing hhy \u\'^. We obtain by the embedding of W 
into L°° 

II T^ I l2ll ^ ^11 T.' I l2l| 

\w 

|2|| , ^111^ ..|2 



|/\ *|mP||loc, ^C||A'*|wpi 



^C||/Ci|Ui|||u|^||v^ + C||/C2|U-|||«rilLi 

2 

L^nw- 



<Ch"2 



So we have 

^\\uo\\l^+C\\u\\1^^,^^T. 
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Moreover, always by applying Lemma 12.21 with h = \u\ we obtain 

mu){t)\\w ^ \\uo\\w + I \\K*\uWT)\\w\\u{T)\\wdT 
Jo 

^ \\uo\\w + C f ||u(T)||i.nH'dT, 

"'0 

and thus 

mu)\\L^w^\\uo\\w + C\\u\\l^^-.^^T. 

Finally we have 

\\'^{u)\\L^L^nw ^ \\uo\\L^nw + C||iio|li2nw^- 
By reducing sufficiently T (depending on ||uo||i2p,^^) we get for all u G X, 
||*&('«)(t)||L«=L2 ^ 2||uo||L2nw- Moreover, for u,v E X we have 

mu) - Hv)\\L^L-nw ^ \\Giu) - G{v)\\l^l2 + ||G(«) - G{v)\\l^w. 
For some positive constant C we have 

\\Giu)-^G{v)h^L^^C\\giu)^giv)\\L.^^. 

< G\\g{u) - g{v)\\L^L^nwT. 
Moreover 

\\Giu) - G{v)\\w ^ f \\giu)ir)-giv){r)\\wdr, 
■Jo 
and 

\\G{u) - G(w)||l~w s$ ll.g(u) - giv)\\L^wT 

^ hiu) - g{v)\\L^L^nwT. 
By replacing /i, /2 by u, v respectively in Lemma [2^ we obtain 
mu) - <^{v)\\L^L2r,w ^ C\\giu) - g{v)\\L^L2r,wT 

Choosing T possibly smaller (still depending on llMolU^nH') we deduce that $ is a 
contraction from X to X. Thus <& has a unique fixed point u E X and Theorem 
O follows. D 

3. Derivation of the approximate solution 



We consider the rescaled version of (|2.1I) : 

=-2 

(3.1) iedtu" + —Au" = eX{K * \u''\^)u'. 
We seek an approximation of solutions to p.ip in the form 

<pp(i,2:)^^a,(t,a;)e^*^(*'-)/^ 

We begin by proceeding formally. We plug the ansatz above into p.ip . This yields 

2 2 

(3.2) tedtulpj, + y A<pp - s{K * K^pPX^^ = ^ e'^Z^ + eiW + r^), 

fc=0 



with 
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and 



(3.3) Z| = i^Aa 






Other terms appear due to the presence of the nonhncarity which are 
(3.4) W = - ^{K * Y, \aif)aje'*'/', 

and 



(3.5) r' = -Y,{K* J2 {akale'^^^-'f"^/'))aje' 



0j/e 






We aim to ehminate all equal powers of e. Hence, by setting Zq = we obtain the 
eikonal equation, whose solution is explicitly given by 

(3.6) 0j(t,a-) ==Kj -a;- -|kj|2. 

Next, we set Zf + W^ = without including r'^ . This latter will constitute the 
first term error, the second will be Zf. We obtain for all j £ N 

(3.7) dtttj + Kj ■ Vflj = —i{K * 2^ k^P)'iji fli(O) = otj, 

where we have used the fact that A^j = 0. 

Lemma 3.1. The transport equation p.7p with initial amplitudes 

(aj)jgN G ii^(IR'') admits a unique global- in-time solutions 

a = (aj)jgN G C{[0, oo[; E{M.)), which can be written in the form 

(3.8) aj(t,x) = aj{x - iKj)e*^^(*''=\ 
where 

(3.9) Sj{t,x)^- (A'*^|a£(x+(r-i)Kj -rK£)p) dr. 

Proof. We multiply p.7p with TiJ. We obtain 

(3.10) aJdtOj + Kj ■ (ajVaj) = -i{K * ^ |a£p)|ajp G iR. 

But 

2Re{ajdtaj + kj • (ajVaj)) = {dt + kj ■ V)|ajp. 
So, we deduce from (|3.10p that (9t + kj ■ V)|ajp = 0, which yields 

\aj{t,x)f = \ajix~tKj)f. 
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and p.Sp follows for some real function Sj. To determine Sj we inject p.8|) into 
dSH). We get 

i{dt + Kj ■ \7)Sj{t, x)aj{x — Ikj) = —{{K * N^ I'^ji^ ~ tne)\'^)ctj{x — tnj)- 

It suffices to impose 

dt{Sj(t,x + tKj)) = -K * ^ |aj(2^ - i(Ki? - Kj))P, 

fSN 

which yields 

Sj{t, X + tKj) = — / K * y. I'^ji^ ~ T{Ke — Kj))p dr, 
and finally wc have 

Sj{t,x) = — / K*y \aj{x + {t ^ t)Kj) ^ TKg\ dr. 

Global in time existence of Uj 's is not trivial from their explicit formula like we 
can see in the following. We rewrite p.7p in its integral form 



(3.11) aj{t,x) ^ aj{x — tKj) + / JV{a)j {t,x + {t — t)Kj)dT, 

Jo 

where the nonlinearity Af is given by 

AA(a), =-z(x*^|a,rta,. 

For a G E{M.'^) we have by definition 

iiA/'(«)iu = Eii-^Kii>' = E E ii5''-A/'KiUw. 

Let 1 77 1 ^ n. Leibnitz formula yields 

||5W(a),|U.ni^=ll9''((if*El«^l')«0lli^nH' 

^cj2 iia^*E^'i«^i')ii^ii«^ii^ 
^ciia,iiy(Eii^iiu^ E ii^'i«^nu^ 

+ II/C2IIL0. ^ iia'^ia.HiH. 

eeN 
^C\\a,\\Y\\a\\l. 
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Finally we have 

\\N-ia)\\E ^ C\\a\\l. 
This shows that Af{a) defines a continuous mapping from E^ to E and by the 
standard Cauchy-Lipschitz theorem for the ordinary differential equations, a 
local-in-time existence results immediately follows. Now, we have to show that 
the solution a(t) = (aj(i))jeN is global in time. Let [0,Tniax[ be the maximal time 
interval where {aj)j^fq is defined. From (|3.11l) we have for JTyj ^ n and t e [0,Tniax[ 

d"aj{t) = d'laji- - tKj) + ( d"J\f{a)j{T, • + (r - t)Kj) dr 

Jo 

ft 



Taking 1? r\W norm yields 

^ Wo.,\W-nw + C(||/Ci|Ui + \\lC2\\L^)Y.\\a,\\l f \\a,{T)\\Y dr 

^W^o^.h^-nw + C f ||a,(r)||ydT, 
Jo 
and so 

hM\Y = E ll^^WIU^nvv ^ \\o^Ay+C f \\a,{T)\\Y dr. 

„i^„ "'0 



By Gronwall lemma 



\\a,my^\\a,\\ye^\ 



for some positive constant C independent of j and t. After summing with respect 
to j we obtain 

\\ait)\\E SC \\a\\Ee^'. 
The growth of ||a(i)||£; is at most exponential, and thus can not explode at finite 
time. We deduce that T„iax = oo and a = {aj)j^fii £ C([0, oo[; E{M.''')). 

D 
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4. Estimations on the remainder 
Wc will estimate in L^ DW the term r"^ and Z|. To this end, we assume 

Proposition 4.1. Let r^ be defined by p.Sp with the plane-wave phases <j)j given 
by (|3.6p . We have the following bound: 



(4.1) VWi^.^w^CS-^-^Mle'''^ 

(4.2) m\\L-nw ^ Me, 
where 6 is like defined in Theorem ] 1. 11 

Proof. Inequality (|4.2|) is obvious. We recall the definition of r^: 

r^ = -"^{K* Y^ afea7e*("^'=-'^^)/^)aje*'^^/^ 

k ^ e 
We estimate W norm of r^, we obtain 



rlw 



Y,iK* J2 afcale*^*''"'*'^^')aje'*^'^'l|i 



jGN k,e em 

k i: I 



<Ell^* E akale'^*'-'^''^''\\w\\aj\\w 



jGN fe,£ e : 



< ^ ||i^*afea?e^(^''-^^)/^||w^||a 



Let us look more closely to the W norm of the convolution above. For h 
sufficiently regular and w G R'', we denote I^{x) = K * (6e"^'^/'^). Remark that 

T^ = {2Tif/'^K be^'^-^/^ 

e 
So, if we take the W norm of I^ we obtain 

wnw = iifiiLi 

= C||A-(«--) (i--) H:--)U, 

+ c\\ic,{^--) (e--) b{---)\\,. 

^ C||/Ci|Ui II (^ - -)'" 5(. - -)|U^, sup U-- 

+ cii(e-^)>-^)iu.sup(ier-(e-^^-" 
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By Peetre inequality we have for all m < d — 7 and |C| > 1 






eI (0 \ £ 



and thus 



Moreover 



UJ 



and 






.-^>%f..(=) .c. 






So, for m = d — "f 

Hence, for allkj€N,k^£ 

\\K * ame^"^ \\w ^ CII/CilUi || (0''"'^ ^A\l^ l^fc " «,|'^-^e'^-^ 
+ C\\1C2\\l^ II {i)"-'' ^IIli \ku - K.^'-^e"-^ 
^ C^^-'^||/Ci|Ui|| (O'^^^IU^^e'^-^ 

Remark that 

II (^} ""'afcOfllL^ =^ II (0"afcaZ||L» 

<c ^ ||a^Ka^)IU- 

<C 5^ ||9''(afeal)|Ui 

?y I ^ n k/ 1 ^ ^^ 

^ C'llafcllylla^llF- 

Moreover 

II (0 "''afcojllii ^ II (^"afea^llii 

^C^ ||9^^K^)|Ui 
s:c^ ||9''Ka?)|k 
<C5] WakWw J2 \\9^<^i\\w 

|?7|^n |77|^n 

^C||afc||r||a,||y. 
Thus wc obtain 
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\\K*akirie^^^\\w^C5'<-\\\lCi\\L^+C\\K.2\\L^)\\ak\\Y\\aA\Ye''-'<. 
In view of the above inequality, we deduce 

VWw^ J2 C6''-''mi\\L^+C\\IC2\\L^^)\\a,\\Y\\a,\\Ye'-''Y.\\a,\\Y 

fcSN feN jGN 

< C<5''-'^(||/Ci|Ui +C||/C2||L..)||a|||£'*-^. 
To control L^ norm of r^ it suffices to remark that 



JSN 


k.l €N 
k^l 






^ y, \\K H 


- > ; flfeoZe^^^' 


.-4>i)h- 


' U°= "i L^ 


jeN 








^c>; 


\\K * akole''-'^''- 


>f)/el 


k) 'llajlly 


k,e et* 






jSN 



and in the same way as previously we get 



Finally 



WWL^nw ^ CP'\\K^\W + ||/C2||L^)||a|||£^-^ 



5. Justification of the approach 



D 



In this section wc show that the solution u| is a good approximation of the exact 
solution in the sense of Theorem 11.11 We assume that assumptions of Theorem 
11.11 are verified. At this stage we have constructed an approximate solution 

in C{[0,oo[;E{M.'^)) where the corresponding 0j's are like constructed in p.6p and 
the profiles are given by Lemma 13.11 Let T^ > the time of local existence of the 
exact solution u^ to p.ip given by the fixed point argument. We consider a fixed 
T > and we introduce the error between the exact solution u"^ and the 
approximate solution u'^pp : 

w'^ =u^ ~ ulpp. 

We have in particular 

ulpp<-C{%T],L^nW), 
so there exists R > independent of e such that 
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Since w^ € C([0, T"^], L^ n W) and w^(0) = 0, there exists i"^ such that 
(5.1) \\w'it)\\L2nw^R, 

for t e [0,t"]. So long as (|5.ip holds, wc infer 

Ik'WIIi^nw s; |A| / II {g{u' + w') - g{ulpp)) Ml^^nw dr 

"'0 

+ |A|£ / \\ZI{T)\\L2nwdT+ |A| / ||i?(T)|U.nH'dT 
Jo Jo 

ft rt 



^C f ||u;-'(T)||i2ni^dT + C£ f dT + Ce'^"'' f dr 
Jo Jo Jo 

^C f \\w'ir)\\L2nwdT + Ce^ f dr, 
Jo Jo 



where we have used Lemma 12.21 and Proposition 231 with /3 = min{l, d — 7}. 
Gronwall lemma implies that so long as (|5.ip holds, 

lk'(i)||L^niv<C(e^*-l)e^ 

for all t G [0,i^]. Choosing e g]0,£o] with £0 sufficiently small, we see that (j5.ip 
remains true for t € [0,T] and Theorem 11.11 follows . 
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